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ABSTRACT 

With the successful launch of the Swift Gamma-ray Burst Explorer, a rich trove of early X-ray 
afterglow data has been collected by its on-board X-Ray Telescope (XRT). Some interesting features 
are emerging, including a distinct rapidly decaying component preceding the conventional afterglow 
component in many sources, a shallow decay component before the more "normal" decay component 
observed in a good fraction of GRBs (e.g. GRB 050128, GRB 050315, GRB 050319, and GRB 050401), 
and X-ray flares in nearly half of the afterglows (e.g. GRB 050406, GRB G5G5G2B, GRB 050607, 
and GRB 050724). These interesting early afterglow signatures reveal valuable and unprecedented 
information about GRBs, including the prompt emission - afterglow transition, GRB emission site, 
central engine activity, forward-reverse shock physics, and the GRB immediate environment. In this 
paper, we systematically analyze the possible physical processes that shape the properties of the early 
X-ray afterglow lightcurves, and use the data to constrain various models. We suggest that the steep 
decay component is consistent with the tail emission of the prompt gamma-ray bursts and/or of the 
X-ray flares. This provides strong evidence that the prompt emission and afterglow emission are 
likely two distinct components, supporting the internal origin of the GRB prompt emission. The 
shallow decay segment observed in a group of GRBs suggests that very likely the forward shock keeps 
being refreshed for some time. This might be caused either by a long-lived central engine, or by 
a wide distribution of the shell Lorentz factors, or else possibly by the deceleration of a Poynting 
flux dominated flow. X-ray flares suggest that the GRB central engine is very likely still active after 
the prompt gamma-ray emission is over, but with a reduced activity at later times. In some cases, 
the central engine activity even extends to days after the burst triggers. Analyses of early X-ray 
afterglow data reveal that GRBs are indeed highly relativistic events, and that early afterglow data of 
many bursts, starting from the beginning of the XRT observations, are consistent with the afterglow 
emission from an interstellar medium (ISM) environment. 

Subject headings: gamma rays: bursts - X-rays: theory - shock waves - radiation mechanisms: non- 
thermal 



1. INTRODUCTION 

With the successful launch of the Swift Gamma-ray 
Burst Explorer, an era of systematic, multi-wavelength 
observations of GRB early afterglows has been ushered 
in. Very early optical/IR detections have been made with 
ground-based telescopes (Akerlof et al. 1999; Fox et al. 
2003; Li et al. 2003; Blake et al. 2005; Vestrand et al. 
2005) before and during the initial operation of Swift. In 
the X-ray band, some evidence of the early afterglows has 
been collected earlier (e.g. Piro et al. 1998; Giblin et al. 
1998; Burenin et al. 1999; Piro et al. 2005). However, 
it is the operation of the Swift X-Ray Telescope (XRT) 
that offers the possibility to unveil the final gap between 
the prompt emission and the late afterglow stage. 

There has been widespread expectation that the early 
X-ray observations could answer a series of core questions 
in GRB studies. What is the connection between the 
GRB prompt emission and the afterglow? Are prompt 
emission and afterglow both from the external shock 



(Meszaros & Rees 1993; Dermer & Mitman 1999) or do 
they come from different locations [i.e. prompt emis- 
sion from the internal shocks (Rees & Meszaros 1994; 
Paczynski & Xu 1994), while the afterglow comes from 
the external shock (Meszaros & Rees 1997a; Sari et al. 
1998)]? Does the central engine become dormant after 
the burst is over? What is the immediate environment of 
the burst, an ISM or a wind? Are there density clumps 
surrounding the GRB progenitor? What is the role of 
the reverse shock? What is the initial Lorentz factor of 
the fireball? 

All these questions can be at least partially answered 
with the early X-ray afterglow data, sometimes in com- 
bination with the prompt gamma-ray data and the early 
optical/IR afterglow data. Although early afterglow 
lightcurves have been extensively modeled in the optical 
band (mainly driven by the observations and by the the- 
oretical argument that the reverse shock emission com- 
ponent plays an important role in the optical band, e.g. 
Meszaros & Rees 1997a; Sari & Piran 1999; Meszaros & 
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Rees 1999; Kobayashi 2000; Kobayashi & Zhang 2003a,b; 
Zhang et al. 2003; Wei 2003; Wu et al. 2003; Li et 
al. 2003b; Fan et al. 2004a; Zhang & Kobayashi 2005; 
Fan et al. 2005a,b; Nakar & Piran 2004; McMahon et 
al. 2004), possible early X-ray afterglow signatures have 
been only sparsely studied (e.g. Kumar & Panaitescu 
2000a; Kobayashi et al. 2005a; Fan & Wei 2005). 

In its first six months of operations, the Swift X-Ray 
Telescope (XRT) has already accumulated a rich col- 
lection of early afterglow features in the X-ray band. 
The XRT is a sensitive broad-band (0.2-10 keV) imager, 
which can be promptly slewed to GRB targets triggered 
by the Burst Alert Telescope (BAT) within 1-2 minutes 
(Burrows et al. 2005a). It is therefore an ideal instru- 
ment to study the transition between the GRB prompt 
emission and the very early X-ray afterglow. The follow- 
ing features are all detected by XRT in a good sample 
of bursts (see, e.g. Chincarini et al. 2005; Nousek et al. 
2005 for a collection of data), reflecting some common 
underlying physics of GRBs. 

1. In most cases (e.g. GRB 050126 and GRB 
050219A), a steep decay is observed up to several hun- 
dred seconds after the burst trigger, which is followed 
by a more conventional shallower decay (Tagliaferri et 
al. 2005; Goad et al. 2005). This conclusion is drawn 
by choosing the GRB trigger time as the zero time 
point {to)- At later times (e.g. {t — to) ^ Tgo, where 
Tgo is the duration of the GRB), the afterglow decay 
slope dlnF^/d\n{t — to) is essentially independent on 
the adopted to- However, at early times (e.g. (t — to) not 
much larger than Tgo, the measured decay slope could be 
very sensitive to the assumed ^o- Tagliaferri et al. (2005) 
explored the to effect and concluded that the two distinct 
lightcurve segments are likely intrinsic rather than due 
to a poor choice of to- Furthermore, in some cases, the 
steep decay segment also has a different spectral index 
(e.g. for GRB 050319, Cusumano et al. 2005). Usu- 
ally it also connects to the spectral extrapolation of the 
BAT prompt emission lightcurve smoothly (Barthelmy 
et al. 2005a). All these facts indicate that the steep de- 
cay component is a distinct new component that may be 
unrelated to the conventional afterglow component. 

2. In a good fraction of GRBs (e.g. GRB 050128, 
Campanaet al. 2005; GRB 050315, Vaughan et al. 2005; 
GRB 050319, Cusumano et al. 2005; and GRB 050401, 
de Pasquale et al. 2005), the early X-ray afterglow is 
characterized by a shallow to "normal" transition. Dur- 
ing the transition the spectral index is not changed. The 
decay slope after the break (e.g. ~ —1-2) is consistent 
with the standard ISM afterglow model, while the decay 
slope before the break is usually much shallower (e.g. 
- -0.5). 

3. In some cases (e.g. GRB 050315, Vaughan et al. 
2005), a further steepening is detected after the "normal" 
decay phase, which is consistent with a jet break. 

4. Bright X-ray flares have been detected in the early 
X-ray lightcurves of nearly one half of the burst popu- 
lation (e.g. GRB 050406, GRB 050202B, Burrows et al. 
2005; Romano et al. 2005; Falcone et al. 2005). In par- 
ticular, the X-ray afterglow of the short-hard burst GRB 
050724 also shows at least three flares (Barthelmy et al. 
2005b). The flares typically happen hundreds of seconds 
after the trigger or earlier, but in some cases they occur 
around a day (e.g. GRB 050502B, Falcone et al. 2005; 



GRB 050724, Barthelmy et al. 2005b). The amphtudes 
of the flares are usually larger than the underlying after- 
glow component by a factor of several (e.g. a factor of 
6 in GRB 050406, Burrows et al. 2005; Romano et al. 
2005), but can be much larger (e.g. ^ 500 in the case of 
GRB 050202B, Burrows et al. 2005; Falcone et al. 2005). 
A similar feature was evident for GRB 011121 detected 
by BeppoSAX (Piro et al. 2005). 

In summarizing the current X-ray afterglow data, one 
can tentatively draw a synthetic cartoon lightcurve in the 
X-ray band, which consists of 5 components (see Figure 
nj: I. an initial steep decay (with a typical slope ~ — 3 or 
steeper); II. a shallower-than-normal decay (with a typ- 
ical slope ^ —0.5); III. a normal decay (with a typical 
slope ^ —1-2); IV. a late steeper decay (with a typi- 
cal slope ~ —2); and V. one or more X-ray flares. We 
note that Nousek et al (2005) also arrived at the similar 
schematic diagram that includes the segments I, II and 
III in our cartoon picture (see their Fig. 3). Limited by 
the quality of the data, the current analyses indicate that 
the spectral indices remain unchanged in the segments II, 
III and IV, with a typical value of /3x ~ 1 {Fx oc v"^^ , 
Nousek et al. 2005). In some bursts, the segments I 
and II have different spectral indices (e.g. GRB 050319, 
Cusumano et al. 2005). In some cases, a time-evolution 
of the spectral index has been detected (e.g. in the giant 
flare of GRB 050502B, Falcone et al. 2005). In Fig.l, 
we have indicated the typical temporal index for each 
segment. Throughout the paper the transition times 
between the adjacent segments for the four lightcurve 
segments are denoted as tfci, tf,2, and if,3, respectively. 
In this paper, we systematically study the physical pro- 
cesses that shape an early X-ray lightcurve, and discuss 
possible theoretical interpretations to the above phenom- 
ena. In 321 we discuss the GRB tail emission arising 
from high angular latitude relative to the viewing direc- 
tion, which takes longer to reach the observer due to 
the extra distance it travels, as a conical (or spherical) 
shell suddenly stops shining. This is the so-called "curva- 
ture effect" (Fenimore et al. 1996; Kumar & Panaitescu 
2000a; Dermer 2004; Dyks et al. 2005). In we review 
the main emission properties from the external forward 
shock region in the X-ray band, summarizing the tempo- 
ral and spectral indices expected in the X-ray band for 
both the ISM and the wind models. Furthermore, we also 
discuss the case of a continuously refreshed shock as well 
as its three possible physical mechanisms. Several case 
studies are investigated to reveal an intriguing refreshed 
shock phase commonly present in many bursts. In f^l we 
briefly discuss whether and how the reverse shock emis- 
sion would influence the X-ray band emission. In f^l we 
explore various mechanisms that might give rise to the 
X-ray flares observed in many bursts (e.g. GRB 050406 
and GRB 050202b), and conclude that the phenomenon 
is best interpreted as due to a late central engine activity. 
Our conclusions are summarized in iJHJ 

2. GRB TAIL EMISSION AND THE CURVATURE EFFECT 

2.1. GRB tail emission 

The temporal bridge between the GRB prompt emis- 
sion and the afterglow emission is essential for revealing 
whether the prompt emission and the afterglow originate 
from the same component. The earliest GRB relativis- 
tic shock model invoked the external shock as the site for 
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Fig. 1. — A synthetic cartoon X-ray lightcurve based on the 
observational data from the Swift XRT. The phase "0" denotes 
the prompt emission. Four power law lightcurve segments together 
with a flaring component are identified in the afterglow phase. The 
segments I and III are most common, and they are marked with 
solid lines. Other three components are only observed in a fraction 
of bursts, so they are marked as dashed lines. Typical temporal 
indices in the four segments are indicated in the figure. The spec- 
tral indices remain unchanged for segments II, III and IV, with a 
typical value of fix ~ 1 (Fx oc The segment I sometimes 

has a softer spectrum (e.g. fix ~ 1.5), but in some other cases 
it has a similar spectral index as the other three segments. The 
flares (segment V) have similar spectra as the segment I, and time 
evolution of the spectral index during the flares has been observed 
in some bursts (e.g. GRB 050502B). 



prompt 7-ray emission (Rees & Meszaros 1992; Meszaros 
& Rccs 1993). The rapid variability observed in many 
GRBs is in great contrast with the intuitive expectations 
in the external shock model, which generally predicts a 
smooth burst temporal profile, and it has been argued 
that the radiative efficiency is too low for the model so 
that a much larger energy budget is required (Sari & Pi- 
ran 1997). Dermer & Mitman (1999, 2003) argued that if 
the GRB ambient medium is sufficiently clumpy, an ex- 
ternal shock GRB model could reproduce the observed 
variability with a high energy efficiency. Within such a 
picture, the prompt emission and the afterglow originate 
from the same component (i.e. the external shock), and 
it is expected that the two emission components are likely 
smoothly connected in the early phase. 

On the other hand, it is now commonly believed that 
GRB prompt emission originates from some "internal 
processes", i.e. the 7-rays are emitted before the fire- 
ball is decelerated by the ambient medium. The most 
widely discussed model is the internal shock model (Rees 
& Meszaros 1994; Paczynski & Xu 1994; Kobayashi et al. 
1997; Daigne & Mochkovitch 1998; Fan et al. 2004b). 
Alternatively, the internal emission may be caused by 
dissipation of strong magnetic fields (e.g. Drenkhahn 
& Spruit 2002), or Comptonization of the photospheric 
emission (Rees & Meszaros 2005). Within such scenar- 
ios, there exist two distinct temporal episodes dominated 
by the prompt emission and the afterglow, respectively, 
since the latter is emitted at a much larger distance from 



the central engine when the fireball is decelerated. Gen- 
erally one should expect a flux contrast between these 
two episodes. 

Before the Swift era, no solid observation was available 
to finally differentiate both scenarios, and evidence in 
favor of each scenario had been collected (see e.g. Zhang 
& Meszaros 2004 for a review). It is one of the major 
tasks of Swift to pin down the emission site of the GRB 
prompt emission. 

If the prompt emission and the afterglow arise from dif- 
ferent emission sites, as is expected in the internal shock 
(or similar) scenario, and if the prompt emission fiux level 
is much higher than the afterglow emission fiux level, one 
expects to sec a steeply decaying lightcurve during the 
transition from the prompt emission phase to the after- 
glow phase. Such a steep decay is due to the so-called 
"curvature effect" (e.g. Kumar & Panaitescu 2000a; Der- 
mer 2004; Dyks et al. 2005; Panaitescu et al. 2005). In 
principle, such an effect also applies to the tail emission 
of the X-ray flares (Burrows et al. 2005). Hereafter we 
generally define such an emission component as "GRB 
tail emission" . 

2.2. Curvature effect 

2.2.1. The simplest case 

For a conical jet with an opening angle 0j, emission 
from the same radius Rcr but from different viewing lat- 
itudes 9 {6 < 9j) would reach the observer at differ- 
ent times. Even if the emission is turned off instan- 
taneously, due to the propagation effect the observer 
would receive the emitted photons at the angle 9 at 
t = {1 + z){Rcr/c){9'^/2). Such a tail emission thus lasts 
for a duration of 



tta^l = (1 + z){Rcr/c){9f/2) ~ 330 S 



if the line of sight is not too close to the jet edge. 

We consider the simplest case of a jet moving with a 
constant bulk Lorentz factor T (or a constant velocity 
?,'). The electrons arc shock- heated up to a radius Rcr, 
beyond which no fresh shocked electrons are injected, 
and the already heated electrons cool rapidly. The co- 
moving emission frequency v' is boosted to = Vi^' in 
the observer's frame, where V = [r(l — u cos 6'/c)]~^ is 
the Doppler factor, which is 7? 2r for 9 ^ l/F, and 
V - 2/(r6'2) for6' > l/F. Since t oc 9'^, one gets D oc t'^ 
for 6* > l/F. 

The observed fiux Fi, is related to the comoving surface 
brightness L^, by 




-2-/3 



(2) 



where f3 is the observed spectral index around the ob- 
served frequency v, and the last proportionality is valid 
for l/F 9 < 9j. With the standard convention 
Fi, oc v~l^t~°', one has the well-known result for the cur- 
vature effect (e.g. Kumar & Panaitescu 2000a; Dermer 
2004; Fan & Wei 2005; Dyks et al. 2005; Panaitescu et 
al. 2005) 

a = 2 + /3 . (3) 
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Fig. 2. — The curvature effect for a decelerating ejecta. In 
tlie calculation, the radiation from the fluid elements on the equal 
arrival time surface for R > Rcut are cut out, and the radiation 
from all other fluid elements is integrated. Here Rcut = 10^^,3 X 
cm are the radii at which the radiation is assumed 
to terminate abruptly. Notice the insensitivity of the decay slope 
(a = 2 + /3) on the cut-oflf radius. The following parameters are 
adopted: the isotropic kinetic energy Eiso = 10^^ ergs, the initial 
Lorentz factor Fq = 240, the jet half-opening angle 6j = 0.1, the 
ambient density n = Icm"'', the redshift z = 1 (the luminosity 
distance = 2.2 X lO^^cm), the electron spectral index p = 2.3, 
and the electron and magnetic equipartition parameters = 0.1 
and €b = 0.01, respectively. 



2.2.2. Emission from a decelerating fireball 

In reality, the Lorentz factor of the shell could be de- 
creasing right before the sudden cessation of the emis- 
sion. This is valid for the external shock case in the de- 
celeration phase, or even in the internal shock case. We 
perform numerical calculations to investigate such an ef- 
fect^. The curvature effect for the sudden switching-off 
of radiation in a decelerating outflow is presented Fig- 
ure El in which we plot the X-ray lightcurve from an 
expanding jet blastwave with Eiso = lO^^ergs, Tq = 240, 
9j = 0.1, n = lcm~'^, and z — 1. We manually turn 
off the radiation at radii {Rcut) of 10^^ cm, 3 x 10^^ cm, 
and 6 x 10^^ cm, respectively, and investigate the subse- 

^ In this paper, we use two codes to perform numerical calcula- 
tions. The first code was developed by J. Dyks (Dyks et al. 2005). 
The code can deal with the afterglow emission of an outflow with 
an arbitrary axisymmetric structure and an arbitrary observer's 
viewing direction. The dynamics of the radial outflow is modeled 
similar to Granot & Kumar (2003). Only synchrotron radiation is 
taken into account at the moment and synchrotron self-absorption 
(which is irrelevant for X-ray emission) is ignored. The synchrotron 
emissivity in the comoving frame is calculated by integration of the 
synchrotron spectrum over the electron energy distribution. The 
latter is calculated by solving the continuity equation with the 
power-law source function Q = K-y^'P , normalized by a local injec- 
tion rate (Moderski et al. 2000). The electrons cool down through 
synchrotron radiation and adiabatic expansion. All kinematic ef- 
fects which affect the observed flux (e.g. Doppler boost, propaga- 
tion time effects) have been taken into account rigorously, following 
Salmonson (2003). The second code was developed by Y. Z. Fan. 
This is an afterglow code developed from the dynamical model of 
Huang et al. (2000), and has been used in several previous studies 
(e.g. Fan et al. 2004a; Fan et al. 2005a). The latest addition is to 
also include the kinetic evolution of the electron distribution (Mod- 
erski et al. 2000). We have used both codes in various calculations 
in this paper, and the results are consistent with each other. The 
figures presented in this paper are all generated from the second 
code. 



quent curvature effect. The radiation from fluid elements 
at i? < Rcut has been rigorously integrated. The result 
indicates that the a = 2 + (3 conclusion is essentially 
unchanged. 

2.2.3. Jet structure effects 

Another interesting issue is the jet structure. In prin- 
ciple, GRB jets could be structured (Zhang & Meszaros 
2002a; Rossi et al. 2002; Zhang et al. 2004a; Kumar & 
Granot 2003). Since the curvature effect allows one to see 
the high-latitude emission directly, an interesting ques- 
tion is whether the decay slope associated with the cur- 
vature effect depends on the unknown jet structure. We 
have investigated this effect with the first code, and find 
that for a relativistic outflow the temporal slope of the 
curvature effect is largely insensitive to the jet structure 
as long as the viewing angle is not far off the bright beam. 
The main reason is that the decrease of flux because of 
the curvature effect occurs on a much shorter timescale 
than that for the jet structure to take effect. For a spec- 
tral index — /3 (/^ cx the flux decreases by m orders 
of magnitude after a time of icrv = lO^^'^^'^^Hcr, where 
tcr is the observer time at which the curvature effect be- 
gan. For a typical /3 1, the flux drops by one order of 
magnitude after a short time tcrv ^ 2tcr- A drop of three 
orders of magnitude occurs in no more than a decade in 
time. On the other hand, the observer can perceive the 
switch-off of emissivity at an angle 6 measured from the 
line of sight at a time tg ~ [\ + {9TY]tcr- One can see that 
the structure of the outflow must have a typical angular 
scale smaller than S/F in order to affect the observed flux 
before lO^cr- For F > 10^, the parameters of the outflow 
would have to vary strongly on a scale smaller than one 
degree. Nonetheless, the effect of the jet struture would 
start to play a noticeable role if the line of sight is out- 
side the bright beam. Detailed calculations are presented 
elsewhere (Dyks et al. 2005). 

2.2.4. Factors leading to deviations from the a — 2 + f3 

relation 

Almost all the Swift XRT early afterglow lightcurves 
are categorized by a steep decay component followed by a 
more "normal" decaying afterglow lightcurve (Tagliaferri 
et al. 2005; Nousek et al. 2005) - see Segment 1 in 
FigCJ In most of these cases, the measured a and /3 
values in this rapidly decaying component are close to the 
a = 2 + (3 relation, but in most cases, they do not match 
completely. This does not invalidate the curvature effect 
interpretation, however, since in principle the following 
factors would lead to deviations from the simple a = 2+/3 
law. 

1. The time zero point (to) effect. In GRB stud- 
ies, the afterglow lightcurves are plotted in the log-log 
scale, with to = defined as the trigger time of the 
burst. When discussing the late afterglows, shifting to 
by the order of the burst duration Tgp does not make 
much difference. When discussing the early afterglow 
and its connection to the prompt emission, however, the 
decay power law index (d In F,y/d In t) is very sensitive to 
the to value one chooses. Correctly choosing to is there- 
fore essential to derive the correct temporal decay index 
a (see e.g. Kobayashi et al. 2005b for a detailed dis- 
cussion about the to issue). In the case of the internal 
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1 (a) 

Fig. 3. — The effect of to on the lightcurves. The same internal 
shock 7-ray pulse is calculated, but is assigned to three different 
ejection times t^j. The dotted, solid, and dash-dotted lines are for 
tcj = 0.0, 1.0, 10.0s, respectively. The following parameters are 
adopted to calculate the internal shock pulses: the pulse luminosity 
ipulsG = lO^^ergs s~^, the variability time scale St = 0.1s, 9j = 
0.2rad, to = 0.5, = 0.1, p = 2.5, and z = I. 



shock model, the case is straightforward. The observed 
pulses essentially track the behavior of the central engine 
(Kobayashi et al. 1997). Each pulse marks another re- 
start of the central engine, so that to should be re-defined 
for each pulse when the curvature effect of that pulse is 
considered. Keeping the same to as the beginning of the 
first pulse (i.e. the GRB trigger) would inevitably lead to 
a false, very steep power law decay for later pulses. Fig- 
ure Ogives an example to show this point. When a series 
of shells are successively ejected, each pulse will be fol- 
lowed by its tail emission due to the curvature effect, but 
most of these tails are buried under the main emission of 
the next pulse. The observed curvature effect is only the 
tail emission of the very last pulse. As a result, properly 
shifting ta is essential to interpret the steep decay com- 
ponent observed in XRT bursts. 2. The superposition 
effect. The observed steep-to-shallow transition in the 
early phase of XRT bursts (Tagliaferri et al. 2005) sug- 
gests that by the end of the tail emission, the fireball is 
already decelerated, and the forward shock emission also 
contributes to the X-ray band. As a result, the observed 
steep decay should also include the contribution from the 
forward shock. Assuming the later has a temporal decay 
index ~w, the X-ray flux at the early phase should read 

where A and B are constants, and tQ i, tg ^ are the time 
zero points for the steep decay component (presumably 
of the internal origin) and for the shallow decay com- 
ponent (presumably of the external origin), respectively. 
In the intermediate regime between the two power-law 
segments, both components are important, and the ob- 
served a should be shallower than 2 + P during the steep 
decaying phase. This effect flattens the decay instead of 
steepening it. 

3. If with the above two adjustments, the observed a is 
still steeper than 2 + P, one can draw the conclusion that 
the solid angle of the emitting region is comparable to 
or smaller than l/F. This would correspond to a patchy 



shell (Kumar & Piran 2000a) or a mini-jet (Yamazaki et 
al. 2004). A caveat on such an interpretation is that 
the probability for the line of sight sitting right on top 
of such a very narrow patch/mini-jet is very small. As a 
result, this model can not interpret an effect that seems 
to be a general property of X-ray afterglows. 

4. If with the first two adjustments, the observed a 
is flatter than 2 + f3 but is still much steeper than that 
expected from a forward shock model, there could be 
two possibilities. One is that the emission is still from 
the internal dissipation of energy but the emission in the 
observational band does not cease abruptly. This is rel- 
evant when the observational band is below the cooling 
frequency. The adiabatic cooling therefore gives a decay 
slope of (1-1- 3/3/2) rather than ~ {2 + (3) (e.g. Sari 
& Piran 1999; Zhang et al. 2003). The second possibil- 
ity is that one is looking at a structured jet (Zhang & 
Meszaros 2002a; Rossi et al. 2002), with the line of sight 
significantly off-axis. The curvature effect in such a con- 
figuration typically gives a flatter decay slope than 2 + (3 
(Dyks et al. 2005). This is particularly relevant for X- 
ray rich GRBs or X-ray flashes for which a large viewing 
angle is usually expected (Zhang et al. 2004; Yamazaki 
et al. 2004) . Further analyses of XRT data suggest that 
at least in some GRBs, the decay slope is shallower than 
2 + f3 (O'Brien et al. 2005). The above two possibilities 
are in particular relevant for these bursts. 

We suggest that most of the rapid-decay lightcurves 
observed by the Swift XRT may be interpreted as GRB 
(or X-ray flare) tail emission through the curvature ef- 
fect, with the first two adjustments discussed above. In 
order to test this hypothesis, after the submission of this 
paper we have performed more detailed data analyses on 
a large sample of XRT bursts (Liang et al. 2005). By 
assuming that the decay slope should be 2 + (3, we search 
the appropriate to that allows such an assumption to be 
satisfied. It is found that to is usually at the beginning of 
the last pulse (for the steep decay following the prompt 
emission) or at the beginning of the X-ray flare (for steep 
decay following flares). This fact strongly suggests that 
the curvature effect is likely to be the correct interpreta- 
tion, and at the same time lends strong support to the 
internal origin of the prompt emission and X-ray flares 
(see §§2.3 & 5.7 for the arguments in favor of the internal 
models for both components). 

Another potential test of the curvature effect is to 
search for a correlation between the spectral peak en- 
ergy {Epk) and the flux at the peak (Fpk) (Dermer 2004). 
This requires a well measured Epk in the XRT band. In 
most cases, the XRT spectrum is consistent with a sin- 
gle power law. More detailed analyses on future bright 
X-ray flares are desirable to perform such a test. 

2.3. Theoretical implications 

The current Swift XRT observations of the early rapid- 
to-shallow decay transition of the X-ray lightcurves 
(Tagliaferri et al. 2005), when interpreted as the cur- 
vature effect, have profound implications for the under- 
standing of the GRB phenomenon. 

1. It gives a direct observational proof that the GRB 
prompt emission very likely comes from a different site 
than the afterglow emission. This suggests that the 
emission either comes from the internal shocks due to 
the collisions among many engine-ejected shells (Rees & 
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Meszaros 1994; Paczynski & Xu 1994), or is due to mag- 
netic or other dissipation processes at a radius smaller 
than the fireball deceleration radius (e.g. Drenkhahn & 
Spruit 2002; Rees & Meszaros 2005). In both scenar- 
ios, the energy dissipation region is well inside the region 
where the deceleration of the whole fireball occurs. 

2. An interesting fact is that in most cases, after the 
prompt emission, the X-ray emission level (that spec- 
trally extrapolated from the BAT data) drops by sev- 
eral orders of magnitude (through the curvature effect, 
in our interpretation) before "landing" on the after- 
glow emission level. One could roughly estimate the 
expected "drop-off". The flux level in the XRT band 
during the prompt phase could be roughly estimated as 
^prompt ^ (E^^,^jT,o){ExKT/Epf+^, whcrc is 
the isotropic energy of the emitted gamma-rays, Tgo is 
the duration of the burst, Exkt ~ 5 keV is the typical 
energy in the XRT band, Ep 100 keV is the typical 
peak energy in the GRB spectrum, and ^ — 1 is the 
low-energy spectral index for a Band-spectrum (Band et 
al. 1993). Assuming that the X-ray band for the after- 
glow emission is above both the typical synchrotron fre- 
quency Vm and the cooling frequency Vc (which is usually 
the case for the ISM model, see eqs.[7], jS|), the X-ray af- 
terglow flux level can be estimated as (e.g. Freedman 
& Waxman 2001) F^^^ cx €eEiso/t, where Eiso is the 
isotropic energy of the afterglow kinetic energy, and Ce is 
the electron equipartition parameter in the shock^. The 
flux contrast can be estimated as 



j-iprompt 



E~j 



Ei, 



t 



{ExYa/Ep) 



a+2 



(5) 



For typical parameters, one has {Exkt I Ep)°'^'^ ^ 0.05 
and ee ~ 0.1, so that the term in the bracket < 1. Al- 
though t > Tgo would generally suggest that pP^°™P*' 
should be higher than -F"^, the large contrast between 
the two components observed in many bursts is usually 
not accounted for unless E^^i^o is (much) larger than 
Eiso- This refers to a very high apparent GRB radi- 
ation efficiency - even higher than the one estimated 
using the late X-ray afterglow data (Lloyd-Ronning & 
Zhang 2004)^^. The commonly-invoked internal shock 
model predicts a low emission efficiency (e.g. Panaitescu 
et al. 1999; Kumar 1999). Understanding such a high 
apparent radiation efficiency is therefore desirable (see 
e.g. Beloborodov 2000; Kobayashi & Sari 2001). 

3. The common stcep-to-shallow transition feature in- 
dicates that the fireball has already been decelerated at 
the time when the GRB tail emission fades. Otherwise, 
one would see an initially rising lightcurve peaking at the 
fireball deceleration time. This fact alone sets a lower 
limit to the intial Lorentz factor of the fireball, since the 
deceleration time tdec must be earlier than the transition 

^ When the synchrotron self-Compton process dominates the 
cooUng, the discussion could be more compHcated. 

^ This could be attributed to the shallow decay injection phase 
(segment II in Fig^as discussed in §3.2. Because of the injection, 
the effective Ei^o in the early epochs is smaller than that in the 
later epochs. As a result, a larger pP^^^P^ _ _F°^ contrast is 
expected. 



TABLE f 

Constraints on the initial Lorentz factors of several 
GRBs. 



GRB 


z 


ibi(s) 




To 




050126.. 


1.290'' 


~ 110 


0.77 


> 120(»?T.n) 


-1/8 


050315.. 


1.949=: 


~ 400 


2.77 


> 100{?7T.n) 


-1/8 


050319.. 


3.240^ 


~ 400 


5.12 


> 120{»7T.n) 


-1/8 


050401.. 


2.900== 


~ 130 


27.49 


> 220(»?Tn) 


-1/8 



^'Chincarini et al. 2005 
''Berger et al. 2005 

Kelson & Berger 2005 
'^Fynbo et al. 2005a 
<=Fynbo et al. 2005b 



time t\,i. The numerical expression is 



''61,2 



1 + 2 



3/8 



(6) 



where E^^iso is the isotropic gamma-ray energy (which 
is an observable if the redshift z is known), 77-y = 
Ej^iso/Eiso is a conversion factor between the isotropic 
afterglow energy Eiso and Ejjso- Throughout the paper, 
the convention Qx = is adopted in cgs units. Ap- 

plying the method to the bursts with measured z (Chin- 
carini et al. 2005), we get the lower limits og Fq for 
several bursts (Table Given the weak dependence on 
the unknown parameters (i.e. (yy^n)"^/^), we conclude 
that the data suggest that GRBs are highly relativis- 
tic events with typical Lorentz factors higher than 100. 
This is an independent method, as compared with previ- 
ous ones using the high energy spectrum (e.g. Baring & 
Harding 1997; Lithwick & Sari 2001), the reverse shock 
data (Sari & Piran 1999; Wang et al. 2000; Zhang et 
al. 2003), and the superluminal expansion of the radio 
afterglow source image (Waxman et al. 1998). 

3. FORWARD SHOCK EMISSION 

After the rapid fading of the GRB tail emission, usually 
the forward shock emission component gives the main 
contribution to the early X-ray afterglow lightcurves. 
The lightcurve shape depends on the density profile of 
the ambient medium (i.e. ISM or wind). In the "stan- 
dard" case (i.e. adiabatic evolution with promp injec- 
tion of energy), the fireball energy is essentially constant 
during the deceleration phase. The bulk Lorentz factor 
F cx foj. the ISM case and F oc i?-^/^ foj. ^^le wind 

case. When the bulk Lorentz factor F is larger than , 
where 9j is the jet opening angle (or the viewing angle of 
a structured jet), the system is simply determined by the 
ratio of the isotropic afterglow energy Eiso and the am- 
bient density n (or the A parameter in the wind model). 
Such a "normal" decay phase corresponds to the segment 
III in the synthetic lightcurve (Fig. When F becomes 
smaller than 9~^, the lightcurve steepens because of the 
combination of the jet edge effect and the possible side- 
ways expansion effect (Rhoads 1999; Sari et al. 1999; 
Panaitescu & Meszaros 1999). The bulk Lorentz factor 
decreases exponentially with radius. This is known as a 
"jet break" , and the post-break segment corresponds to 
the segment IV in Fig. ^ 

During the early evolution of the fireball, the forward 
shock may be continuously refreshed with additional en- 



7 



ergy. This could be either because of a continuous op- 
eration of the central engine (Dai & Lu 1998a; Zhang & 
Meszaros 2001; Dai 2004), or because of a power law dis- 
tribution of the Lorentz factors in the ejecta that results 
in slower ejecta catching up with the decelerated fire- 
ball at later times (Rees & Meszaros 1998; Panaitescu et 
al. 1998; Kumar & Piran 2000b; Sari & Meszaros 2000; 
Zhang & Meszaros 2002b), or because of the transfer- 
ing of the Poynting-flux energy to the medium when a 
Poynting-flux dominated flow is decelerated (e.g. Zhang 
& Kobayashi 2005). The cannonical XRT lightcurve 
(Fignj indeed shows a shallow decay phase (segment II), 
which we argue is due to continuous energy injection. 

3.1. Standard afterglow models 

For the convenience of the later discussion, we summa- 
rize the "standard" early forward shock X-ray afterglow 
properties as follows. 

The ISM model (e.g. Sari et al. 1998). The typi- 
cal synchrotron frequency and the cooling frequency are 
= 6.5x1014 Hz e]^^_^el_iElih~^^^[il+z)/2]^/^, and 
1^, = 2.5 X IQie Hz (1 + y)-2e^3/2_E-V2^-i^-i/2j(^ ^ 

z)/2]~i/^, respectively, where E is the isotropic kinetic 
energy of the fireball, n is the ISM density, Ce and es are 
shock equipartition parameters for electrons and mag- 
netic fields, respectively, Y is the energy ratio between 
the inverse Compton component and the synchrotron 
component, z is the redshift, t is the observer's time. 
Both frequencies decrease with time. The time interval 
for and to cross the XRT energy band (0.5-10 keV) 
from above can be expressed as 

t„, = (4 - 30) s e]i'_,et!l,Eli' (7) 

te = (0.1-40) s {l + Y)-%%E^,'n-' (^i±^) '(8) 

The epoch when the fireball switches from fast cooling 
(i^c < Vrn) to slow cooling (vc > Vm) IS defined by requir- 
ing Um — i^c, which reads 

W = 26s (l + r)2e| _2<-i£^52ri(^i^^ . (9) 

For comparison, the time when the fireball is decelerated 
(thin shell case) is given by 

= 180s(i?52/n)V3ro-f , (10) 

where Fq is the initial Lorentz factor of the fireball. We 
can see that for typical parameters, the XRT band is 
already in the regime of vx > max ( I'm, Vc) when deceler- 
ation starts. Also, the blast wave evolution has usually 
entered the slow cooling regime where the radiative losses 
are not important^. Under such conditions, the temporal 

In certain parameter regimes, the condition t^ec < tmc could 
be satisfied, and in the temporal regime t < tmc, the blastwave is 
in the fast cooling (uc < Vm) domain, so the radiative loss could 
afli'ect the blastwave dynamics (e.g. Bottcher & Dermer 2000; Wu 
et al. 2005). 



index is^ ax — (3p — 2)/4 ^ 1.15, and the spectral in- 
dex is Px = p/2 ~ 1.1 (photon index 2.1). The relation 
between the two indices is ax = {^I3x — l)/2. In view 
that tc is very sensitive to e^, one can still get the regime 
Vm < vx < Vc if e_B is large (say, around 0.1). In such a 
case, ax = 3(p - l)/4 ~ 0.9, f3x = {p - l)/2 ~ 0.6, and 
ax = (3/2)/3x- The above two cases have been com- 
monly observed in early X-ray afterglows of many Swift 
bursts (e.g. segment HI in Fig. suggesting that the 
fireball shock model can successfully interpret the general 
properties of GRB afterglows, and that most GRBs are 
born in a constant density medium. This conclusion is 
consistent with previous analyses (Panaitescu & Kumar 
2002; Yost et al. 2003). 

The wind model (e.g. Chevalier & Li 2000). The typi- 
cal synchrotron frequency and the cooling frequency are 

,yra = 1.3 X 1015 Hz el/,i2<-i^^52\'^'[(l + ^)/2]'/' 
(notice that in the wind case has the same pa- 
rameter dependences as the ISM case but is larger by 
a factor of ^ 2, e.g. Dai & Lu 1998b), and Vc = 

6.8x1012 Hz il + Y)-\-%lEli'A-Hl^'[il + z)/2]-^/^ 
where = {Mw / 4.TrVw) / {5 x 10" g cm"i) is the typi- 
cal wind parameter, Mw is the mass loss rate, and Vw is 
the wind velocity. A distinguishing property of the wind 
model is that Vc increases with time. Similarly, the time 
interval for i/m to cross the XRT band from above is 

U = (6.6 - 49) s eZ^CAi' (^) , dD 

and the time interval for Vc to cross the band from below 
is 



tc = (3.1 X 10" - 1.2 X 10") s (1 + Y)'^el _2Er^^At 




= (0.3-120)s(l + y)V5,_2i?52'<-3(^) • ( 

The critical time for the fast-slow cooling transition is 

tmc = 1.4 X lO-i s (1 + y)ee,-ieB,_2A 

= 14s (l + r)ee,-ieB,-2A,_3 . (13) 

The deceleration time is essentially the duration of the 
burst, i.e. tdec ~ T = 20 s T1.3, since for typical pa- 
rameters, the wind model is the so-called "thick shell" 
case (e.g. Chevalier & Li 2000; Kobayashi & Zhang 
2003b). For a typical wind parameter A^ ^ 1, the 
X-ray lightcurve is very simple. For t < i„j, one has 
Vc < vx < so that ax — 1/4, (3x = 1/2 and 
ax = (1 — /3x)/2. When t > tm, during most of the ob- 
servational time of interest, one has vx > max(i^m, Vc), 
so that ax = (3p-2)/2- 1.15, /3jf = p/2 = 1.1 (photon 
index 2.1), and ax = (3/3x — l)/2. The switching be- 
tween the fast cooling and slow cooling regimes does not 
influence the temporal and spectral indices in the X-ray 
band. Only when A^, < 0.01, i.e. tc falls into the range 
of observational interest, does a new temporal/spectral 
domain appear. When t > tc, one has Vm < vx < Vc, 

^ Here and below, the convention F^{X) oc t~°'^ u~l^^ is 
adopted since in the X-ray band both the temporal and the spectral 
indices are negative. Also p = 2.2 is adopted for typical numerical 
values. 
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ax = (3p- l)/4 - 1.4, Px = (p- l)/2 ~ 0.6 (photon in- 
dex 1.6), and ax = (3/?x + l)/2. Such a feature has been 
used to interpret GRB 050128 (Campana et al. 2005)^. 
If j4» is not much smaller than unity, the blastwave is in 
the fast cooling regime, and radiative losses could be sub- 
stantial (Bottcher & Dermer 2000). A detailed analysis 
has been presented in Wu et al. (2005). 

The jet model (e.g. Rhoads 1999; Sari et al. 1999). Af- 
ter the jet break, the temporal decay index is predicted to 
be ax = P- This is derived by assuming significant side- 
ways expansion. This result is independent on whether 
the X-ray band is below or above Uc, and whether the 
medium is an ISM or a stellar wind. For the latter, the 
time scale for the lightcurve to achieve the asymptotic 
—p index is typically longer than that in the ISM case 
(e.g. Kumar & Panaitescu 2000b; Gou et al. 2001). 

All the above discussions apply for the case of p > 2. 
For p < 2, the case could be different. Dai & Cheng 
(2001) proposed one scenario to deal with the case of 
p < 2, while Panaitescu & Kumar (2002) extended the 
treatment of p > 2 case to the p < 2 regime. 

3.2. Refreshed shock models 

If there is significant continuous energy injection into 
the fireball during the deceleration phase, the forward 
shock keeps being "refreshed" , so that it decelerates less 
rapidly than in the standard case. The bulk Lorentz 
factor of the fireball decays more slowly than F oc 
R^'^/^{R^^/^) for the ISM case (the wind case), respec- 
tively. 

There are three possible physical origins for the re- 
freshed shocks. 

1. The central engine itself is longer lasting, e.g. be- 
having as 

m = Lo{t/U)-^ . (14) 

The dynamical evolution and the radiation signature of 
such a system has been discussed in detail in Zhang & 
Meszaros (2001). A specific model for such an injection 
case, i.e. the energy injection from the initial spin down 
from a millesecond pulsar (preferably a millisecond mag- 
netar) was discussed in that paper and earlier in Dai & 
Lu (1998a). In such a specific model, g = is required 
according to the spin-down law. Alternatively, the con- 
tinued engine activity could be due to continued infall 
onto a central black hole, resulting in the time depen- 
dence eo. (11411 ^. In general, for an adiabatic fireball, the 
injection would modify the blastwave dynamics as long 
as q < 1 (Zhang & Meszaros 2001). The energy in the 
fireball increases with time as Eiso oc t^~'^, so that 

r oc R^'^VT^ (xt-^,R (xt^, ISM (15) 
F oc cx oc t^,wind (16) 

It is then straightforward to work out the temporal in- 
dices for various temporal regimes. 

The ISM model. The typical synchrotron frequency 
Vrn oc F^7ei? OC F^ OC t^(2+g)/2^ ^j^g synchrotron cooling 

As discussed in il3.2l after collecting more data, we now believe 
that the shallow-to-normal decay observed in GRB 050128 is more 
likely due to the transition from the energy injection phase to the 
standard phase (without injection). 

The black hole - torus system typically has q = 5/3 at later 
times (MacFadyen et al. 2001; Janiuk et al. 2004), which has no 
effect on the blastwave evolution. 



frequency oc T-^B-^t-^ oc T-H-^ oc t(9-2)/2^ 
the peak flux density Fi^^„iax oc N^BT oc i^^'^, where 
i? oc F is the comoving magnetic field strength, 7e oc 
F is the typical electron Lorentz factor in the shocked 
region, and oc R^ is the total number of the emitting 
electrons. The temporal indices a for various spectral 
regimes and their relationships with the spectral indices 
a{(3) are listed in Tabled 

The wind model. In the wind case, the ambient density 
is n oc R~^, where R is the radial distance of the shock 
front to the central source. The typical synchrotron fre- 
quency Vm oc F^7e-B oc F'^B oc i^(2+(?)/2^ ^Yie synchrotron 
cooling frequency i^c oc T^^B^-^t^^ oc T^'^t^^ oc i(2-g)/2^ 
and the peak fiux density F^^max oc N^BT oc F^ oc 
where B oc TR~^ is the comoving magnetic field 
strength, and Ne oc Ris the total number of the emitting 
electrons. The temporal indices a for various spectral 
regimes and their relationships with the spectral indices 
a{P) are listed in TableEl 

In order for the central engine to continuously feed 
the blast wave, the Lorentz factor of the continuous flow 
must be (much) larger than that of the blast wave. It 
could be a Poynting-flux-dominated flow. This is not 
difficult to satisfy since the blast wave keeps decelerat- 
ing. There could be a reverse shock propagating into the 
continuous ejecta, but the radiation signature of the re- 
verse shock is typically not in the X-ray band (e.g. Zhang 
& Meszaros 2001). 

2. The central engine activity may be brief (e.g. as 
brief as the prompt emission itself, but at the end of the 
prompt phase, the ejecta has a range of Lorentz factors, 
e.g., the amount of ejected mass moving with Lorentz fac- 
tors greater than 7 is (Rees & Meszaros 1998; Panaitescu 
et al. 1998; Sari & Meszaros 2000) 



M(> 7) oc 7 



(17) 



The total energy in the fireball increases as Eiso oc 



7I-" oc Pi-^ so that 



F oc R-— oc t-^,Roc t^,ISM 



(18) 



Foci? !+'■ oc t 3+= , i? oc , wind (19) 

One can then work out the temporal decay indices in 
various spectral regimes (e.g. Rees & Meszaros 1998; 
Sari & Meszaros 2000). Alternatively, for each s value, 
one can find an effective q value that mimics the s effect, 
or vice versa. This gives 



10- 7? 

2 + q ' 
4-3(7 



10 - 2s 

4 



ISM 



. wind 



(20) 



(21) 



q 3 

In Table [3 the explicit s-dependences are not listed, but 
they could be inferred from eas. (|20|l and (|21() . 

In this second scenario, the central engine need not last 
long. All the material could be ejected promptly. The 
continuous injection is due to the different velocities of 
the ejecta. Initially as the blast wave moves with high 
speed, the slower ejecta lag behind and have no effect on 
the blastwave evlolution. They later progressively pile up 
onto the blast wave as the latter decelerates. Only when 
s > 1 does one expect a change in the fireball dynamics. 
This corresponds to g < 1. For q = 0.5, one gets s = 2.6 
for the ISM case and s 5 for the wind case. 
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TABLE 2 

Temporal index a and spectral index (3 in various 

AFTERGLOW MODELS. 





/9 


no injection 
a 


a(/3) 


injection 
a 








ISM 


slow cooling 
















!^ < f m 
Um < I' < 
V > l/c 


_ 1 

(0.65) 
f (1-15) 


_ 1 

3(P-1) ^^ n\ 


3/3 

3/3 — 1 
"= 2 


(2p-6) + (p+3)9 


(0.3) 
(0.7) 


a={q- 

a= (q- 

q — 2 


1) + ^-'\'>>' 
1) + '^+''^'' 


4 

3p-2 
4 


(1.2) 


4 

{2p-4) + (p+2)<j 
4 


1 (2 + 9)^ 

9. 


ISM 


fast cooling 
















U < Uc 

Uc < U < Urn. 

U > Vm 


1 

3 

1 

2 

f (1.15) 


1 
6 

1 

4 

3p-2 


(1.2) 


" 2 
" 2 


22=8 (_o 8) 
^ (-0.1) 

(2p— 4)-|-(p-|-2)(j 


(0.7) 


" = (<?- 

a = (q- 
^ 2 


1) + (^7^^ 

1) + (^7>'^ 
1 (2-1-9)^ 


4 




4 




2 


Wind 


slow cooling 
















I'm < ^ < 
V > Vc 


_ 1 

£=i (0.65) 
f (1-15) 




3p-l 


(1.5) 
(1.2) 


2 

3/3-1 

a- 


2^ (-0.2) 

{2p-2) + (p+l)9 


(1.1) 
(0.7) 


«= f + 

"=f + 

9-2 


(2-H9)/3 

(2+\)/3 


4 

3p-2 
4 


(2p-4) + (p+2)5 
4 


1 h+l)0 

2 


Wind 


fast cooling 
















U < Vc 

Vc < U < I'm 

l'> ^m 


1 

3 

1 (1-15) 


2 
3 
1 
4 

3p-2 
4 


(1.2) 


1-/3 

" 2 

3/3-1 
ry = — ^ 


^'+''' (0.5) 
^ (-0.1) 

(2p-4)-Kp-|-2)g 
4 


(0.7) 


" 2 

" 2 

9-2 
" 9 


(2-9)/3 

(2-9)/3 

1 ?2-|-g)/3 
2 



Note. — This is the extension of the Table 1 of Zhang & 
Meszaros (2004), with the inclusion of the cases of energy injection. 
The case of p < 2 is not included, and the self-absorption effect 
is not discussed. Notice that a different convention Fi, oc t~°'i'~^ 
is adopted here (in comparison to that used in Zhang & Meszaros 
2004), mainly because both the temporal index and the spectral in- 
dex are generally negative in the X-ray band. The temporal indices 
with energy injection are valid only for q < I, and they reduce to 
the standard case (without energy injection, e.g. Sari et al. 1998, 
Chevalier & Li 2000) when q = I. For g > 1 the expressions are no 
longer valid, and the standard model applies. An injection case due 
to pulsar spindown corresponds to </ = (Dai & Lu 1998a; Zhang 
& Meszaros 2001). Recent Swift XRT data are generally consistent 
with q ~ 0.5. The numerical values quoted in parentheses axe for 
p = 2.3 and q = 0.5. 
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3. The energy injection is also brief, but the outflow 
has a significant fraction of Poynting flux (e.g. Usov 
1992; Thompson 1994; Meszaros & Rees 1997b; Lyutikov 
& Blandford 2005). Assigning a parameter cr for the 
outflow, which is the ratio between the Poynting flux and 
baryonic kinetic energy flux, Zhang & Kobayashi (2005) 
modeled the reverse shock emission from ejecta with an 
arbitrary a value. They found that during the crossing of 
the reverse shock, the Poynting energy is not transferred 
to the ambient medium. The Poynting energy (roughly 
by a factor of a) is expected to be transferred to the 
medium (and hence, to the afterglow emission) after the 
reverse shock disappears. Zhang & Kobayashi (2005) 
suggest that the transfer is delayed with respect to the 
traditional case of cr = 0. The energy transfer process, 
however, is poorly studied so that one does not have a 
handy conversion relation with the q value derived in the 
first scenario. 

3.3. Case studies 

In this subsection, we discuss several Swift GRBs with 
well- monitored early afterglow data detected by XRT. 
The notations for the break times and the temporal 
slopes are per those marked on Fig. 

GRB 050128 (Campana et al. 2005): The hghtcurve 
can be fitted by a broken power law with the break 
time at tb2 — 1472j;2goS. The temporal decay indices 

before and after the break are a2 — 0.27^0 :^0 and 
as = l.SO^Q j^g, respectively. The spectral indices be- 
fore and after the break are essentially unchanged, i.e. 
P2 ~ 0.59 ± 0.08, and /Jg = 0.79 ± 0.11. Campana et 
al. (2005) discussed two interpretations. A jet model re- 
quires a very flat electron spectral index, i.e. p ~ 1.3, as 
well as a change of the spectral domain before and after 
the jet break. Alternatively, the data may be accom- 
modated in a wind model, but one has to assume three 
switches of the spectral regimes during the observational 
gap from 400s to about 4000s. So neither explanation is 
completely satisfactory. By comparing the predicted in- 
dices in Table El the observation may be well-interpreted 
within the ISM model with an initial continuous energy 
injection episode. The segment after the break is consis- 
tent with a standard ISM model for < < i^c, with 
p ~ 2.6. The lightcurve before the break, on the other 
hand, is consistent with an injection model with p 2.2 
and q ~ 0.5 in the same spectral regime. The break time 
is naturally related to the cessation of the injection pro- 
cess, and a slight change of electron spectral index (from 
2.2 to 2.6) is required. From the beginning of the obser- 
vation (100s) to tt,2, the total energy is increased by a 
factor of (1472/100)(i-°-5) - 2.8. 

GRB 050315 (Vaughan et al. 2005): After a steep 
decay (ai = 5.2tQ 4) up to tbi — 308s, the hghtcurve 
shows a flat "plateau" with a temporal index of a2 = 
0.06to °3. It then turns to ^ 0.71 ± 0.04 at = 
1.2^03 X lO'^s. Finally there is a third break at t^s = 
2.5^03 X lO^s, after which the temporal decay index is 

a4 = 2.0to;3- ^^is burst displays all four segments 
presented in Fig^ The spectral indices in segments II, 
III and IV are essentially constant, i.e. P2 = 0.73 ± 
0.11, Ps = 0.79 ± 0.13 and P4 = 0.7toi, respectively. 
Segment III is consistent with an ISM model with vx > 



Vc and p — 1.6, since in this model /3 — p/2 = 0.8, 
a = {ip — 2)/4 = 0.7, in perfect agreement with the 
data. The third temporal break ^(,3 is consistent with a 
jet break. According to Dai & Cheng (2001), the post- 
break temporal index for p < 2 is a = (p -|- 6)/4 = 
1.9, which is also consistent with the observed a4. The 
plateau between 4i and th2 is then due to an energy 
injection in the same ISM model {vx > t'c), with p ~ 1.5 
and q ^ 0.35. The total injected energy is increased by 
a factor of (12000/308)(i""-35) - 11. 

GRB 050319 (Cusumano et al. 2005): After a steep 
decay (ai = 5.53 ± 0.67) up to 4i = (384 ± 22)s, the 
lightcurve shows a shallow decay with a temporal index 
of a2 = 0.54 ± 0.04. It steepens to as = 1.14 ± 0.2 
at tb2 = (2.60 ± 0.70) x lO-^s. The spectral indices in 
segment II and III are /Ja = 0.69 ± 0.06 and P3 = 0.8 ± 
0.08, respectively. Again segment III is well consistent 
with an ISM model for < i^x < i^c with p = 2.6, 
which gives 13 ^ {p - l)/2 = 0.8 and a = (3/2)/3 = 
1.2, in excellent agreement with the data. Interpreting 
the segment II (the shallow decay phase) as the energy 
injection phase, for the same ISM model [v„i < vx < 
Vc), one gets p ^ 2.4 and q ^ 0.6. The total injected 
energy is increased by a factor of (26000/384)'^^^"-^^ ^ 
5.4. The UVOT observations are also consistent with 
such a picture (Mason et al. 2005). 

GRB 050401 (de Pasquale et al. 2005): The early X- 
ray lightcurve is consistent with a broken power law, with 
a2 = 0.63 ± 0.02, ag = 1.41 ± 0.1, and tb2 = 4480l^4oS. 
The spectral indices before and after the break are all 
consistent with /32 ~ /Sa = 0.90 ± 0.03. The a - f3 rela- 
tion does not fit into a simple p < 2 jet model. On the 
other hand, the energy injection model gives a natural 
interpretation. After the break, the lightcurve is consis- 
tent with an ISM model for i^m < i^x < with p — 2.8. 
Before the break, it is consistent with the same model 
with q = 0.5. The total injected energy is increased by a 
factor of > (4480/200)(i-°-5) ~ 4.7. 

The injection signature is also inferred in other bursts 
such as GRB 050117 (HiU et al. 2005) and XRF 050416 
(Sakamoto et al. 2005), where similar conclusions could 
be drawn. The injection model is supported by an inde- 
pendent study of Panaitescu et al. (2005). 

3.4. Theoretical implications 

The following conclusions could be drawn from the 
above case studies. 

1. A common feature of the early X-ray afterglow 
lightcurves is a well defined temporal steepening break. 
A crucial observational fact is that there is essentially 
no spectral variation before and after the break. This 
suggests that the break is of hydrodynamic origin rather 
than due to the crossing of some typical frequencies of the 
synchrotron spectrum in the band. It is worth mention- 
ing that a lightcurve transition similar to the transition 
between segments II and III is expected in a radiative 
fireball (e.g. Bottcher & Dermer 2000), see e.g. Figs. 1& 
2 of Wu et al. (2005). However, that transition is due to 
the crossing of v„i in the observational band. One there- 
fore expects a large spectral variation before and after 
the break, which is inconsistent with the data. Another 
straightforward interpretation would be a jet break, but 
there are three reasons against such an interpretation. 
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First, in all the cases, p < 2 has to be assumed. This is 
in stark contrast to the late jet breaks observed in the 
optical band, which typically have p > 2. Furthermore, 
the a — P relation predicted in the jet model is usually not 
satisfied. Second, the post-break a — (3 relation is usually 
satisified in a standard slow cooling ISM model, with the 
X-ray band either below or above the cooling frequency. 
In such a sense, this segment is quite "normal" . Third, in 
some cases (e.g. GRB 050315), another steepening break 
is observed after this normal segment, which is consis- 
tent with the jet break interpretation. Since only one 
break could be attributed to a jet break, the "shallow- 
to-normal" break must be due to something else. 

2. A natural interpretation of the shallow decay phase 
is to attribute it to a continuous energy injection, so that 
the forward shock is "refreshed" . Three possibilities exist 
to account for the refreshed shock effect i ^'i.2l : a long- 
lived central engine with progressively reduced activities, 
an instantaneous injection with a steep power-law distri- 
bution of the shell Lorentz factors, and the deceleration 
of an instantaneously-injected highly magnetized (high- 
er) flow. In terms of afterglow properties, these possi- 
bilities are degenerate (e.g. the connection between q 
and s) and can not be differentiated. In principle, the 
first scenario may give rise to additional observational 
signatures (e.g. Rees & Meszaros 2000; Gao & Wei 2004, 
2005), which may be used to differentiate the model from 
the others. 

3. Two interesting characteristics during the injection 
phase are that the injection process is rather smooth, and 
that the effective q value is around 0.5. This gives inter- 
esting constraints on the possible physical mechanisms. 
(1) For the scenario of a continuously injecting central en- 
gine (Zhang & Meszaros 2001), the central engine lumi- 
nosity must vary with time smoothly. This is in contrast 
to the conventional GRB central engine which injects en- 
ergy erratically to allow the observed rapid variability in 
the lightcurves. This usually requires two different en- 
ergy components, i.e. one "hot" fireball component that 
leads to the prompt emission and a "cold" Poynting fiux 
component that gives to the smooth injection. A natu- 
ral Poynting fiux component is due to the spin-down of 
a new-born millisecond pulsar (Dai & Lu 1998a; Zhang 
k. Meszaros 2001). However, a straightforward predic- 
tion from such a model is g = 0, not consistent with 
q ~ 0.5 inferred from the observations. Modifications to 
the simplest model are needed. Alternatively, the sys- 
tem may be a long-lived black hole torus system with 
a reducing accretion rate. However, at later times the 
long-term central engine power corresponds to q — 5/3 
(MacFadyen et al. 2001; Janiuk et al. 2004), too steep to 
give an interesting injection signature. It is worth men- 
tioning that in the coUapsar simulations (MacFadyen et 
al. 2001), an extended flat injection episode sometimes 
lasts for ~ 1000s, which could potentially interpret the 
short injection phase of some bursts, but is difficult to 
account for some other bursts whose injection phase is 
much longer. (2) For the scenario of a power-law distri- 
bution of Lorentz factors (Rees & Meszaros 1998), one 
should require that a smooth distribution of Lorentz fac- 
tors is produced after the internal shock phase. In the in- 
ternal shock model, slow shells are indeed expected to fol- 
low the fast shells, but they tend to be discrete and give 
rise to bumpy lightcurves (e.g. Kumar & Piran 2000a) 



especially when the contribution from the reverse shock 
is taken into account (Zhang & Meszaros 2002b). It is 
also unclear how an effective q ~ 0.5 is expected. (3) 
Deceleration of a promptly-ejected Poynting-flux dom- 
inated flow (e.g. Zhang & Kobayashi 2005) naturally 
gives a smooth injection signature observed. Above case 
studies indicate that the injected energy is by a factor 
of several to 10. Within such a picture, the unknown a 
value is about several to 10. However, it is unclear how 
long the delay would be and whether one can account for 
the shallow decay with q ^ 0.5 extending for 10** seconds. 
Further more detailed theoretical modeling is needed to 
test this hypothesis. 

4. Any model needs to interpret the sudden cessation 
of the injection at ti,s- This time has different meanings 
within the three scenarios discussed above. (1) Within 
the long-lived central engine model, this is simply the 
epoch when the injection process ceases. In the pulsar- 
injection model, there is a well-defined time for injec- 
tion to become insignificant (Dai & Lu 1998a; Zhang & 
Meszaros 2001), but within a black- hole-torus injection 
model, such a time is not straightforwardly defined. (2) 
In the varying Lorentz factor scenario, this time corre- 
sponds to a cut-off of the Lorentz factor distribution at 
the low end below which the distribution index s is flat- 
ter than 1 so that they are energetically unimportant. 
This lowest Lorentz factor is defined by 
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A successful model must be able to address a well-defined 
Vm in this model. (3) Within the Poynting flux in- 
jection model, a well-defined time cut-off is expected, 
which corresponds to the epoch when all the Poynting 
energy is transferred to the blastwave. If the shallow de- 
cay is indeed due to Poynting energy transfer, the cut-off 
time {tb2) could be roughly defined by the a parameter 
through a {tb2/tdecY^~'^\ where tdec is the conven- 
tional deceleration time defined by Eiso/i^ + f), when 
only a fraction of (1 -|- cr)~^ energy is transfered to the 
ISM (Zhang & Kobayashi 2005). 

5. Although we have not tried hard to rule out a wind- 
model interpretation, the case studies discussed above 
suggest that the early afterglow data are consistent with 
an ISM model for essentially all the bursts. This conclu- 
sion also applies to other well-studied Swift bursts (e.g. 
GRB 050525a, Blustin et al. 2005). This result is in- 
triguing given that long GRBs are associated with the 
death of massive stars, from which a strong wind is ex- 
pected. Previous analyses using late time afterglow data 
(e.g. Panaitescu & Kumar 2002; Yost et al. 2003) have 
also suggested that most afterglow data are consistent 
with an ISM model rather than a wind model. In order 
to accommodate the data, it has been suggested that the 
wind parameter may be small so that at a late enough 
time the blastwave is already propagating in an ISM (e.g. 
Chevafier et al. 2004). The Swift resufts push the ISM 
model to even earlier epochs (essentially right after the 
deceleration) , and indicate the need for a re-investigation 
of the problem. The epoch shortly before the deaths of 
massive stars is not well studied (Woosley et al. 2003). 
One possibility is that the stellar wind ceases some time 
before the star collapses. Careful analyses of early after- 
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glows of a large sample of long GRBs may shed light on 
the final stage of massive star evolution. 



4. REVERSE SHOCK EMISSION 
4.1. Synchrotron emission 

It is generally believed that a short-lived reverse shock 
exists during the intial deceleration of the fireball and 
gives interesting emission signatures in the early after- 
glow phase. Given the same internal energy in both 
the forward-shocked and the reverse-shocked regions, the 
typical synchrotron frequency for the reverse shock emis- 
sion is typically much lower than that in the forward 
shock region, since the ejecta is much denser than the 
medium. While the early forward shock synchrotron 
emission peaks in X-rays at early times, the reverse shock 
synchrotron emission usually peaks in the optical/IR 
band or even lower (e.g. Meszaros & Rees 1997a; Sari & 
Piran 1999; Kobayashi 2000; Zhang et al. 2003; Zhang 
& Kobayashi 2005). This model has been succussful in 
interpreting the early optical emission from GRB 990123 
(Akerlof et al. 1999; Sari & Piran 1999; Meszaros & Rees 
1999), GRB 021211 (Fox et al. 2003; Li et al. 2003; Wei 
2003), and GRB 041219a (Blake et al. 2005; Vestrand et 
al. 2005; Fan et al. 2005b). As a result, it is expected 
that the reverse shock component has a negligible con- 
tribution in the X-ray band. 

In the above argument, it has been assumed that the 
shock parameters (cg, es and p) are the same in both 
shocks. In reality this might not be the case. In par- 
ticular, the GRB outflow may itself carry a dynamically 
important magnetic field component (or Poynting flux). 
This magnetic field would be shock-compressed and am- 
plified, giving a larger effective es (Fan et al. 2004a, b; 
Zhang & Kobayashi 2005). Since the medium is generally 
not magnetized, it is natural to expect different es values 
in both regions, and a parameter TZb = {f-B,rl^B.fY have 
been used in the reverse shock analysis. It has been found 
that TZb is indeed larger than unity for GRB 990123 and 
GRB 021211 (Zhang et al. 2003; Fan et al. 2002; Kumar 
& Panaitescu 2003; Panaitescu & Kumar 2004; MaMa- 
hon et al. 2004). Hereafter the subscript/superscript 
"f" and "r" represent the forward shock and the reverse 
shock, respectively. According to Zhang & Kobayashi 
(2005), the case of GRB 990123 corresponds to the most 
optimized case with ct ^ 1, so that T^s is the largest. 

More generally, eg and p may also vary in both shocks. 
Fan et al. (2002) performed a detailed fit to the GRB 
990123 data and obtained el = A.7ei and ejj = 400e|j. 
A general treatment therefore requires that we introduce 
one more parameter, i.e. TZe = [{p'^ — 2)/(p'' — l)]/[(p^ — 
2)/(p^ ~ l)](£c/£c)- Following the treatment of Zhang et 
al. (2003), we have the following relations in the thin- 
shell regime (see also Fan & Wei 2005)^ 
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° For an arbitrary a, the treatment becomes more complicated. 
The treatment presented here is generally valid for cr < 1. For 
(T > 1, the reverse shock emission starts to be suppressed (Zhang 
& Kobayashi 2005). Since we are investigating the most optimistic 
condition for the reverse shock contribution, in this paper we adopt 
the standard hydrodynamic treatment which is valid for c < 1. 
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where and V are the Compton parameters for the 
forward and the reverse shock emission components, re- 
spectively; is the reverse shock crossing time, which 
is essentially tdec (eq-CHl) for the thin shell case; Fx is 
the bulk Lorentz factor of the outflow at 734, x ~ 
(Fo/Fx -I- Fx/Fo)/2 is the Lorentz factor of the shocked 
ejecta relative to the unshocked one. 

For typical parameters, both and i/^ are below the 
XRT band (comparing eqs.[7], [S] with eq.^Uj). Accord- 
ing to eas. (|23|l and l|24|l . j/J, i/J^j should be also below the 
XRT band. Following the standard synchrotron emission 
model, we then derive the X-ray fiux ratio of the reverse 
shock and the forward shock components at <x : 
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We can see that for T^s = TZe = 1, — F'', and 
p > 2, one has Fl^ity^) ^ Ft,xitx), since 734, x < 1, 
and Fx ~ Fq in the thin shell case. The reverse shock 
contamination in the X-ray band is therefore not impor- 
tant. The situation changes if we allow higher TZe and 
TZb values. Increasing TZe directly increases the reverse- 
to-forward fiux ratio. Although the dependence on TZb 
is only mild when p is close to 2, a higher TZb sup- 
presses the IC process in the reverse shock region rel- 
ative to that in the forward shock region, so that the 
ratio (1 -I- Y^)/{\ + y) also increases. As a result, as 
TZb ^ and TZe ^ I, the reverse shock synchotron 
component would stick out above the forward shock syn- 
chrotron component, and an X-ray bump is likely to 
emerge (see also Fan & Wei 2005). As a numerical ex- 
ample, taking p = 2.3, Fx « Fo/2 50, TLb = 10, 
TZe = 5, and = 306^3, we get Fl^{ty,) / Fl^{t.) « 6. 
This could potentially explain the X-ray flare (by a factor 
of - 6) detected in GRB 050406 (Burrows et al. 2005; 
P. Romano et al. 2005, in preparation). However, a big 
caveat of such a model is that one expects a very bright 
UV/Optical flash due to the large TZb and TZe involved 
- like the case of GRB 990123. Unless this flash is com- 
pletely suppressed by extinction, the non-detection of 
such a flash in the UVOT band for GRB 050406 strongly 
disfavors such an interpretation. 

4.2. Synchrotron self-Compton emission 

The synchrotron photons in the reverse shock region 
will be scattered by the same electrons that produce 
these photons. The characteristic energy of this com- 
ponent is typically in the 7-ray range. However, under 
some conditions, this synchrotron self-Compton (SSC) 
component would also stand out in the X-ray band, giv- 
ing rise to an X-ray bump in the lightcurve. A detailed 
discussion has been presented in Kobayashi et al. (2005), 
which we do not repeat here. The general conclusion is 
that the SSC component could account for an early X-ray 
flare bump by a factor of several under certain optimized 
conditions. An advantage of this model over the reverse 
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shock synchrotron model is that a bright UV-optical flash 
is avoided. However, this model can not account for a 
flare with a very large contrast (e.g. by a factor of 500, 
as seen in GRB 050502B, Burrows et al. 2005). 

5. MECHANISMS TO PRODUCE EARLY X-RAY FLARES 

XRT observations indicate that X-ray flares are com- 
mon features in the early phase of X-ray afterglows (the 
component V in Figure After the report of the first 
two cases of flares in GRB 050406 and GRB 050502B 
(Burrows et al. 2005), later observations indicate that 
nearly half of long GRBs harbor early X-ray flares (e.g. 
O'Brien et al. 2005). More intriguingly, the early X- 
ray afterglow of the latest localized short GRB 050724 
(Barthelmy et al. 2005b) also revealed flares similar to 
those in the long GRBs (e.g. 050502B). The common 
feature of these flares is that the decay indices after the 
flares are typically very steep, with a 5t/t much smaller 
than unity. In some cases (e.g. GRB050724, Barthelmy 
et al. 2005), the post-flare decay slopes are as steep as 
< —7. In this section we discuss various possible models 
to interpret the flares and conclude that the data re- 
quire that the central engine is active and ejecting these 
episodic flares at later times. 

5.1. Emission from the reverse shock? 

As discussed in 21 synchrotron or SSC emission from 
the reverse shock region could dominate that from the 
forward shock emission in the X-ray band under cer- 
tain conditions. Because of the lack of strong UV-optical 
flares in the UVOT observations, we tentatively rule out 
the reverse shock synchrotron emission model. The pre- 
diction of the lightcurve in the SSC model could poten- 
tially interpret the X-ray flare seen in GRB 050406 (Bur- 
rows et al. 2005), but the predicted amplitude is too low 
to interpret the case of GRB 050502B (Burrows et al. 
2005). In some bursts (e.g. GRB 050607), more than 
one flare are seen in a burst. Although one of these flares 
may be still interpreted as the reverse shock SSC emis- 
sion, an elegant model should interpret these flares by a 
model with the same underlying physics. We tentatively 
conclude that the reverse shock model cannot account 
for most of the X-ray flares detected by the Swift XRT 
in a unified manner. 

5.2. Density clouds surrounding the progenitor? 

Long-duration GRBs are believed to be associated with 
the deaths of massive stars, such as Wolf-Rayet stars. 
According to Woosley et al. (2003), there are no observa- 
tions to constrain the mass loss rate of a Wolf-Rayet star 
during the post-helium phase (100-1000 years before the 
explosion). No stability analyses have been carried out 
to assess whether such stars are stable. At the high end 
of a reasonable range of the mass loss rate, dense clouds 
surrounding a GRB progenitor may exist. For a wind ve- 
locity ^ 100km s~^, the density clumps could occur 
at a radius ^3x10^^ - 3x10 cm. These density clumps 
have been invoked by Dermer & Mitman (1999, 2003) to 
interpret the GRB prompt emission variabilities. 

One immediate question is whether these density 
clumps, if they exist, could give rise to the X-ray flares 
detected by XRT. In order to check this possibility, we 
investigate the following toy model. For simplicity, we 
assume a dense clump extends from i? ^ 3 x lO^^cm 
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Fig. 4. — The X-ray flare powered by a relativistic fireball 
interacting with a dense cloud. The cloud is located at i? = 
(3.0,3.3) X lO^S cm with a density n^iaud = 100 cm-^. The back- 
ground ISM density is n = 1 cm~^. Other parameters include: 
Eiso = 1052ergs, z = 1, Fq = 240, te = 0.1, cb = 0.01, p = 2.3 
and dj =0.1. The dashed line is the X-ray emission contributed 
by the electrons shocked at _R < 3.3 X lO^^cm, and the dotted 
line is the X-ray emission contributed by the electrons shocked at 
-R > 3.3 X lO^^cm. The solid line is the sum of these two compo- 
nents. 



to (i? -I- A) ^ 3.3 X lO^^cm with a number density 
ndoud ~ lOOcm^-^. The background ISM density is taken 
as n ~ 1 cm^"^. Other parameters in the calculation in- 
clude: Ei^o = 10^2gj.gg^ z = 1, Fo = 240, eo = 0.1, 
eB ~ 0.01, p ~ 2.3 and 9-^ — 0.1. No sideways expan- 
sion of the jet is included, which is consistent with pre- 
vious simulations (Kumar & Granot 2003; Cannizzo et 
al. 2004). The X-ray lightcurve is shown in Fig01 the 
general feature of which is reproduced in both codes used 
in this paper. Although the rising phase could be very 
sharp, the decaying slope is rather flat. This is because 
of two effects. First, after entering the dense cloud, the 
blastwave Lorentz factor falls off rapidly. The observed 
time t ~ R/2T^ is hence significantly stretched. The 
dashed line in Fig^represents the observed emission be- 
fore the fireball exits the cloud. Second, after the fireball 
exits the cloud (dotted line in0J, the fireball does not de- 
celerate immediately since the Lorentz factor is already 
too low to be further decelerated in a medium with a low 
density. It is decelerated again when enough medium is 
swept up at a larger radius. The interaction of a fireball 
with density bumps has been studied previously by many 
authors (e.g. Lazzati et al. 2002; Dai & Wu 2003). Our 
detailed calculations suggest that the variation caused by 
density inhomogeneities is generally not very significant 
(see also Ramirez- Ruiz et al. 2005). The lightcurves for 
a blastwave surfing on a density wave should be gener- 
ally quite smooth (cf. Lazzati et al. 2002). The one- 
dimensional model presented here effectively calculates 
the feature when the fireball encounters a high density 
shell. A more realistic model for density clumps should 
include the size of the clump (A). The lightcurve decay 
slope after the peak therefore depends on the comparison 
between F"^ and the angle the clump extends from the 
central engine, A/R. If F^i < A/R, the above calcula- 
tion is still valid since the observer would not notice the 
edge of the clump. The resulting decay slope is rather 
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shallow after the peak, which is in distinct contrast to 
the X-ray flare data. This is the case for the parame- 
ters we use in Fig. 0] For smaller clumps, the lightcurve 
will steepen as F"^ becomes comparable to A/i?. Our 
numerical calculations indeed show such a feature. How- 
ever, the flux does not fall rapidly right after the peak. 
The peak time is typically shortly after the blastwave 
enters the density clump when the Lorentz factor is still 
very high. The epoch when F"^ — A/R is satisfied hap- 
pens much later. The lightcurve therefore still show a 
shallow decay segment before the steep decay compo- 
nent. This is in contrast with the X-ray flare data that 
show rapid decays right after the flare peak (e.g. Bur- 
rows et al. 2005; Falcone et al. 2005; Romano et al. 
2005; O'Brien et al. 2005). One could make the clouds 
small enough, so that F^^ > A/i? is satisfied from the 
very beginning. However, in such a case it is very hard 
to achieve a large contrast between the flares and the un- 
derlying afterglow level (loka et al. 2005). The 500 time 
contrast observed in GRB 050502B is in any case very 
difficult to achieve within all the density clump models 
we have explored. We therefore tentatively conclude that 
the X-ray flares commonly detected in the early X-ray af- 
terglow lightcurves are likely not caused by the putative 
density clouds surrounding the GRB progenitors. The 
lack of these well-modeled features in the data also sug- 
gest that the lumpiness of circumburst medium is, if any, 
rather mild. We note, however, Dermer (2005, in prepa- 
ration) suggests that more detailed 3-D modeling of the 
density clump problem could potentially reproduce the 
observed X-ray flares. 

5.3. Two-component jet? 

In the coUapsar progenitor models, it is expected that 
the ejecta generally has two components, one ultra- 
relativistic component powering the GRB, and another 
moderately relativistic cocoon component (e.g. Zhang, 
et al. 2004b; Meszaros & Rees 2001; Ramirez-Ruiz et 
al. 2002). This gives a physical motivation to the phe- 
nomenological two-component jet model (e.g. Lipunov 
et al. 2001; Berger et al. 2003; Huang et al. 2004). 
An interesting possibility would be whether the X-ray 
flare following the prompt gamma-ray emission (which 
arises from the central relativistic component) is caused 
by the deceleration of the wider mildly-relativistic co- 
coon component as it interacts with the ambient medium. 
A straightforward conclusion is that the wide, off-beam 
component must contain more energy than the narrow 
component in order to give noticeable features in the 
lightcurve. Also the decay after the lightcurve peak of the 
second component should follow the standard afterglow 
model, and the variability time scale satisfies St/t ^ 1. 
The optical lightcurves of the two component jet model 
have been calculated by Huang et al. (2004) and Gra- 
not (2005). FiglSlshows a sample X-ray lightcurve in the 
two-component jet model. It is obvious that this model 
cannot interpret the rapid fall-off observed in the XRT 
X-ray flares, and is therefore ruled out. 

5.4. Patchy jets? 

A related model considers a jet with large energy fluc- 
tuations in the angular direction, so that its energy dis- 
tribution is patchy (Kumar & Piran 2000b). This is a 




Fig. 5. — The X-ray lightcurve powered by a two-component jet. 
The angular range < 9 < 0.1 is the central narrow component, 
which has -Biso,n = lO^^ergs and Fo^n = 240. The wide jet compo- 
nent covers a range of 0.1 < 8 < 0.3, with i?iao,w = 5 X lO^'^ergs 
and Fo,m = 50. The ISM density is n = lcm~^. Other parameters: 
Co = 0.1, EB = 0.01, p = 2.3, and 2 = 1. The line of sight is at 
= 0. The lightcurve peak of the second component corresponds 
to the epoch when the off-beam wide component is decelerated, so 
that its 1/Fiu beam enters the field of view. 




Fig. 6. — The X-ray lightcurve powered by a patchy jet. For 
simplicity, an annular patchy jet is simulated. Following param- 
eters are adopted. Patch 1 (the on-beam jet): < 9 < 0.02, 
Eiso.i = lO^^gj-gs^ ro,i = 240; Patch 2: 0.02 < 9 < 0.04, 
Aso,2 = 5 X lO^^ergs, ro,2 = 50; Patch 3: 0.04 < 9 < 0.06, 
Asoia = lO^^ergs, ro,3 = 240; Patch 4: 0.06 < 9 < 0.08, 
Eiso,4 = 5 X lO^^ergs, Fq = 50; Patch 5: 0.08 < 9 < 0.10, 
Eiso.5 = lO^^ergs and Fo,5 = 240. Other parameters are the same 
as Fig 13 



variation of the two-component jet, and could be approx- 
imated as a multi-component jet. When the 1/F cone of 
different patches enter the field of view, the observed 
lightcurve may present interesting signatures. However, 
the general feature of the two-component jet still applies: 
Only when a patch has a substantial energy compared 
with the on-beam jet would it give a bump feature on 
the lightcurve. After each bump, the afterglow level is 
boosted and would not resume the previous level. The 
variability time scale is also typically St/t ~ 1. Figure 
gives an example. Apparently this model can not ac- 
count for the observed X-ray flares. 
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5.5. Post energy injection into the blastwave? 

In the internal shock model, it is expected that after 
the collisions are finished, the shells are distributed such 
that a shell with a higher Lorentz factor always leads 
a shell with a lower Lorentz factor. As the fast moving 
shell (blastwave) is decelerated, the trailing slow shell will 
catch up with it and inject energy into the blastwave (Ku- 
mar & Piran 2000a). Such an injection also happens if 
the central engine further ejects high-F shells that catch 
up with the decelerating blastwave (Zhang & Meszaros 
2002b). Such a colhsion would give rise to a bump sig- 
nature on the lightcurve. A detailed treatment suggests 
that the overall lightcurve should include emission from 
three components: a forward shock propagating into the 
medium, a second forward shock propagating into the 
blastwave, and a reverse shock propagating into the in- 
jected shell. In the X-ray band, the contribution from 
the reverse shock is negligible. The lightcurve generally 
shows a step-like signature, due to the increase of the to- 
tal energy in the blastwave (Zhang & Meszaros 2002b). 
After the peak, the flux level does not resume the pre- 
vious level since more energy has been injected into the 
fireball. The decay slope after the injection peak follows 
the standard afterglow model, and St/t ^ 1 is expected. 
These are also inconsistent with the data of X-ray flares. 

5.6. Neutron signature? 

Another interesting possibility is whether X-ray flashes 
are the signature of the existence of free neutrons in 
the flrebafl. Derishev et al. (1999), Beloborodov (2003) 
and Rossi et al. (2004) suggested that a baryonic fire- 
ball contains free neutrons, the decay of which would 
leave important imprints on the lightcurve. Fan et al. 
(2005a) modeled the process carefully and calculated the 
lightcurves. According to Fan et al. (2005), the neutron- 
feature is rather smooth in a wind model, and it is hard 
to detect. In the ISM case, on the other hand, a bump 
does exist (in all bands). The physical reason is that 
the trailing proton shell catches up with the decelerated 
neutron-decay products. The physical process is analo- 
gous to the post energy injection effect discussed in ii5.5l 
The amplitude of the fiare is modest, at most a factor 
of several. Since the injection model is not favored, this 
possibility is also disfavored. 

5.7. Late central engine activity 

After ruling out various "external-origin" mechanisms, 
we are left only with the possibility that involves the 
re-activation of the central engine. In this interpreta- 
tion, the X-ray flares share essentially the same origin as 
the prompt gamma-ray emission, i.e. they are caused by 
some "internal" energy dissipation processes which occur 
before the ejecta is decelerated by the ambient medium. 
The leading scenario is the "late" internal shock model, 
which suggests that the central engine ejects more energy 
in the form of an unsteady wind with varying Lorentz 
factors at a late time. These discrete shells collide with 
each other and produce the observed emission. Alterna- 
tively, the late injection could be mainly in the form of 
magnetic fields/Poynting flux, and the X-ray flares are 
due to the intermittent dissipation of the magnetic fields, 
likely through magnetic reconnection events. Fan, Zhang 
& Proga (2005c) argued that at least for the fiares fol- 



lowing short GRBs, the process that powers the flares 
has to be magnetic-origin. 

There are two advantages of the "internal" models over 
the "external" models. 

First, re-starting the central engine equivalently re-set 
the time zero point. In this interpretation, the observed 
flare component and the underlying decaying component 
observed at the same observer time t originate from dif- 
ferent physical sites at different central engine times. Let 
us assume that the initial burst lasts Tgo, that the central 
engine re-activates after a time interval At, and that it 
ejects an unsteady wind with a typical variability time 
scale of 6t. At the observer time t = Tgo + At ^ 1000 
s, the underlying decaying component (external after- 
glow) happens at a distance R^x ~ [r(t)]^ct 5 x 10^^ 
cm, where r(t) is the Lorentz factor of the blastwave at 
the time t. The flare, on the other hand, happens at 
a distance of i?i„ ^ F^cJt ^ (10^^ — lO"^'*) cm, where 
Fq is the initial Lorentz factor of the late-time ejecta. 
According to the clocks attached to the central engine, 
the photons from the external afterglow component are 
emitted at tex ~ Rex/c ^ 1.7 x lO^s, while the pho- 
tons from late central engine activity are emitted at 
Un - t + R,nlc - 103 + i?„/c - (103 - 104)s. Be- 
cause of the relativistic effect, these photons reach the 
observer at exactly the same time t, and superpose onto 
the lightcurve detected by the observer. When plotted as 
a single log — log lightcurve with the origin at the burst 
trigger, a very steep apparent decay slope can be pro- 
duced for a large time shift At (FiglSJ. This naturally 
overcomes the St/t > I constraint encountered by the ex- 
ternal models (loka et al. 2005 and references therein). 

Second, invoking a late central engine activity greatly 
eases the required energy budget. In most of the exter- 
nal models, in order to give rise to a significant bump 
on the lightcurve, the total newly-added energy (either 
from the radial direction - late injection case, or from the 
horizontal direction - patchy jets or multiple-component 
jets) must be at least comparable to the energy that de- 
fines the afterglow emission level. This model therefore 
demands a very large energy budget. For the internal 
model, on the other hand, since the lightcurve is a su- 
perposition of two independent physical components, the 
energy budget is greatly reduced, especially if the bump 
happens at later times when the background afterglow 
level is much lower. For example, for an X-ray lightcurve 
with the decay index of -1 following a burst with dura- 
tion 10 s, a significant flare at ~ 10'*s only requires a 
luminosity slightly larger than 10"^ times of that of the 
prompt emission. This model is therefore very "economi- 
cal" in interpreting the very late (> 10'*s) flares detected 
in some bursts (e.g. GRB 050502B, Falcone et al. 2005; 
and GRB 050724, Barthelmy et al. 2005b). 

Can late central engine activity give rise to softer 
bursts (e.g. X-ray flares as compared with the prompt 
gamma-ray or hard X-ray emission)? According to the 
internal shock model, the peak energy Ep of the syn- 
chrotron spectrum satisfies (Zhang & Meszaros 2002c) 

Ep oc TB' cx L^/V^i cx L^/^T-^6t-^ , (27) 

where L is the luminosity, B' cx (L/F^r^)^/^ is the co- 
moving magnetic field strength, F is the typical Lorentz 
factor of the wind (for internal shock collisions, this F 
is for slow shells), and St is the variability time scale. 
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We can sec that a smaller L, a higher F and a larger 
St would be favorable for a softer burst. The observed 
XRT X-ray flares generally have a smaller luminosity L. 
The lightcurve is smoother with a larger 5t (Burrows et 
al. 2005; Barthelmy et al. 2005b). Also, at later times, 
the environment tends to be cleaner so that P could be 
larger. One therefore naturally expects softer flares at 
later times. 

Can the central engine re-start after being quiescent 
for some time (e.g. ^ lO^s, but sometimes as late as 
~ lO^s)? This is an interesting theoretical problem. The 
collapsar model predicts a central engine time scale of 
minutes to hours (MacFadyen et al. 2001). The star 
may also fragment into several pieces during the col- 
lapse (King et al. 2005). The largest piece collapses 
first onto the black hole, powering a prompt gamma-ray 
burst. Other fragments are intially ejected into elliptical 
orbits, but would eventually fall into the black hole after 
some time, powering the late X-ray flares. Fragmenta- 
tion could also happen within the accretion disk itself 
due to gravitational instabilities (Perna et al. 2005). A 
magnetic-dominated accretion flow could also give rise 
to intermittent accretion flows due to interplay between 
the gravity and the magnetic barrier (Proga & Begelman 
2003; Proga & Zhang 2005). Flares could also occur in 
other central engine scenarios (Dai et al. 2005). The 
rich information collected by the Swift XRT suggests that 
we are getting closer to unraveling of the details of the 
bursting mechanisms of GRBs. More detailed studies are 
called for to unveil the mystery of these explosions. 

An interesting fact (see next subsection for case stud- 
ies) is that the duration of the flares are positively corre- 
lated with the time at which the flare occurs. The later 
the flare, the longer the flare duration. A successful cen- 
tral engine model must be able to address such a peculiar 
behavior. Perna et al. (2005) suggest that if the accre- 
tion disk is fragmented into blobs or otherwise has large 
amplitude density fluctuation at large radii from the cen- 
tral engine, the viscous disk evolution would cause more 
spread for blobs further out from the central engine. This 
gives a natural mechanism for the observed correlation. 
Perna et al. (2005) suggest that gravitational instability 
in the outer part of the disk is likely the origin of the 
density inhomogeneity within the disk. 

5.8. Case studies 

In this subsection, we briefly discuss the X-ray flares 
discovered in several GRBs. 

GRB 050406 & GRB G50502B (Burrows et al. 
2005): These were the first two bursts with flares de- 
tected by XRT. For the case of GRB 050406 whose Tqq 
is ^ 5s, an X-ray flare starts at ^ 150s and reaches a peak 
at ~ 230s. The flux rebrightening at the flare peak is by a 
factor of 6. The rising and the decaying indices are ~ 4.9 
and —5.7, respectively. The total energy emitted during 
the flare is about 10% of that emitted in the prompt 
emission. This suggests that the central engine becomes 
active again after ^ 150s, but with a reduced power. For 
the case of GRB 050502B whose Tgo is ~ 17.5s, a giant 
flare starts at ~ 300s and reaches a peak at ^ 740s. The 
flux rebrightening is by a factor of ~ 500. The decay 
index after the peak is —6. The total energy emit- 
ted during the flare is comparable to that emitted in the 
prompt emission. This suggests that the central engine 



re-starts and ejects a substantial amount of energy. The 
duration of the X-ray flare is much longer that the dura- 
tion of the prompt emission (Tgo), so that the luminosity 
of the flare is much lower. In GRB 050502B, there is yet 
another late flare-like event that peaks at ~ 7 x lO^s. 
The post-peak decay index is ~ —3. This is also consis- 
tent with the case of a late-time central engine activity, 
and the decaying slope is consistent with that expected 
from the curvature effect. The total energy emitted in 
this flare is ~ 10% of that emitted in the giant flare at 
~ 740s. 

GRB 050724 (Barthelmy et al. 2005b): This is a 
short, hard burst with Tgo ~ 3s whose host galaxy is 
an elliptical galaxy, similar to the case of the first Swift- 
localized short burst GRB 050509B (Gehrels et al. 2005). 
The XRT lightcurve reveals rich features which are quite 
similar to the case of GRB 050502B. The XRT obser- 
vation starts at ~ 74s after the trigger, and the early 
XRT lightcurve initially shows a steep decay with a slope 
~ —2. This component is connected to the extrapo- 
lated BAT lightcurve that shows a flare-like event around 
(60 — 80)s. This extended flare-like epoch (including 
the -2 decay component) stops at ~ 200s after which 
the lightcurve decays even more rapidly (with an index 
< — 7). A second, less energetic flare peaks at ~ 300s, 
which is followed by another steep (with index < —7) de- 
cay. A third, significant flare starts at ^ 2 x 10'*s, and the 
decay index after the peak is ~ —2.8. This quite similar 
to those the late flare seen in GRB 050502B. The energy 
emitted during the third flare is comparable to that of 
the second one, and both are several times smaller than 
the energy of the first fiare. All the rapid-decay compo- 
nents following the fiares could be potentially interpreted 
as the high-latitude emission, given a proper shift of the 
zero time point. Since the duration of the third flare is 
very long, the to effect does not affect the decay index 
too much. This is why the decay index —2.8 is quite 
"normal", i.e. consistent with the —2 — (3 prediction. 
The earlier steeper decays (e.g. < —7) are all preceded 
by flares with a sharp increasing phase. Shifting to in 
these cases would lead to signiflcant flattening of the de- 
caying index, which could be still consistent with the 
high-latitude emission. All these discussions also apply 
to the case of GRB 050502B''. 

GRB 011121 (Piro et al. 2005). This BeppoSAX 
burst also indicated a flare-like event around ^ 270s. 
Piro et al. (2005) interpreted the X-ray bump as the 
onset of the afterglow phase. In view of the fact that 
X-ray flares are commonly detected in Swift, bursts, it 
is natural to spectulate that the event is also an X-ray 
flare caused by late central engine activity. Fan & Wei 
(2005) have suggested the late central engine activity and 
performed a detailed case study on this event. 

6. CONCLUSIONS 

During the past several months, the Swift XRT has col- 
lected a rich sample of early X-ray afterglow data. This, 
for the first time, allows us to peer at the final temporal 

gap left by the previous observations and to explore many 
interesting questions of GRB physics. In this paper, we 

^ After submitting this paper, detailed data analyses (Liang et 
al. 2005) confirm these speculations. The late flares in both GRBs 
050724 and 050502B are consistent with the hypothesis that they 
are due to late time central engine activity. 
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have systematically investigated various possible physi- 
cal processes that could give interesting contributions to 
the early X-ray afterglow observations. This includes the 
tail emission of the prompt gamma-ray emission, both 
the forward shock and the reverse shock emission com- 
ponents, refreshed shocks, post energy injection, medium 
density clumps near the burst, angular inhomogeneities 
of the fireball, emission component due to the presence 
of free neutrons, as well as emission from late central en- 
gine activity. We discuss how the above processes might 
leave interesting signatures on the early X-ray afterglow 
lightcurves. 

Based on the XRT data collected so far, we summa- 
rize the salient features and suggest a tentative synthetic 
lightcurve for the X-ray afterglow lightcurves. As shown 
in Fig. ^ the synthetic lightcurve includes five compo- 
nents: I. an intial steep decay component; II. a shallow- 
than-normal decay component; III. a "normal" decay 
component; IV. a post-jet break component; and V. X- 
ray flares. The components I and III appear in almost all 
the bursts. Other three components also commonly ap- 
pear in some bursts. Flares have been detected in nearly 
half of the XRT early lightcurves, and the shallow decay 
segment has also been discovered in a good fraction of 
Swift GRBs. We therefore believe that they represent 
some common underlying physics for GRBs. After com- 
paring data with various physical models, we tentatively 
draw the following conclusions. 

1. The rapid decay component (Tagliaferri et al. 2005) 
commonly observed in the very early afterglow phase 
(which usually has a different spectral slope than the 
late shallow decay components) is very likely the tail 
emission of the prompt gamma-ray bursts or of the early 
X-ray flares. Allowing proper shifting of the time zero 
point and considering the contribution of the underly- 
ing forward shock emission, we speculate that essentially 
all the steep decay cases could be understood in terms 
of the "curvature effect" of the high-latitude emission as 
the emission ceases abruptly. More detailed data analy- 
ses (Liang et al. 2005) support such a speculation. 

2. The transition between the prompt emission and the 
afterglow emission appears to be universally represented 
by a rapid decay followed by a shallower decay, indicating 
that the GRB emission site is very likely different from 
the afterglow site, and that the apparent gamma-ray ef- 
ficiency is very high. 

3. In a good fraction of GRBs (e.g. GRBs 050128, 
050315, 050319, 050401, Campana et al. 2005; Vaughan 
et al. 2005; Cusumano et al. 2005; de Pasquale et al. 
2005), a clear temporal break exists in the early X-ray 
lightcurves. There is no obvious spectral index change 
across the break. The temporal decay index before the 
break is very flat, while that after the break is quite 
"normal" , i.e. is consistent with the standard afterglow 
model for a fireball with constant energy expanding into 
an ISM. We suggest that these breaks are likely not "jet 
breaks" . Rather they mark the cessation of an early con- 
tinuous energy injection phase during which the external 
shock is refreshed. We suggest three possible physical 
mechanisms for the refreshed shocks, i.e. a long-lived 
central engine with a decaying luminosity, a power law 
distribution of the shell Lorentz factors before decelera- 
tion begins, and the deceleration of a high-tr flow. Fur- 
ther studies are needed to better understand this phase. 



4. In the "normal" phase, the data for many bursts 
are consistent with an ISM medium rather than a wind 
medium. This has important implications for under- 
standing the massive star progenitors of long GRBs, in- 
cluding their late evolution stage shortly before explo- 
sion. 

5. Given that most of the shallow-to-normal tran- 
sitions are due to the cessation of the refreshed shock 
phase, the cases with a well-identified jet break are not 
very common. Nonetheless, jet breaks are likely identi- 
fied in some bursts, e.g. GRB 050315 (Vaughan et al. 
2005), GRB 050525 A (Blustin et al. 2005), XRF 050416 
(Sakamoto et al. 2005), and some others. 

6. The X-ray flares detected in nearly half of the Swift 
bursts are most likely due to late central engine activity, 
which results in internal shocks (or similar energy dissi- 
pation events) at later times. It seems to us that there 
is no evidence for the existence of density clumps in the 
GRB neighborhood, and there is no support for strong 
angular inhomogeneities (e.g. two-component jet, patchy 
jets) for the GRB fireball. However, their existence is not 
ruled out. 

7. The similar lightcurves for some long GRBs (e.g. 
GRB050502B, Burrows et al. 2005) and some short 
GRBs (e.g. GRB050724, Barthelmy et al. 2005) indi- 
cate that different progenitor systems may share some 
similarities in the central engine. This might be caused 
by the same underlying physics that controls the hyper- 
accreting accretion disk - a common agent in charge both 
types of systems (Perna et al. 2005). 

The Swift XRT is still rapidly accumulating data on 
the early X-ray afterglows. More careful statistical anal- 
yses of these X-ray data, as well as more detailed case 
studies, also including low frequency data collected by 
UVOT and other ground-based telescopes, would greatly 
imporve our knowledge about GRB prompt emission and 
afterglows, advancing the quest for the final answers to 
the core questions in the study of GRBs. 
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